On the Stability of Linear Discrete-Time Fuzzy Systems 

Gabriele Oliva*, Stefano Panzieri* and Roberto Setola^ 

* University Roma Tre of Rome, Italy, oliva@dia.uniroma3.it, panzieri@uniroma3.it 
t University Campus Bio-Medico of Rome, Italy, r.setola@unicampus.it 



Abstract 

In this paper the linear and stationary Discrete-time systems with state vari- 
ables and dynamic coefficients represented by fuzzy numbers are studied, 
providing some stability criteria, and characterizing the bounds of the set of 
solutions in the case of positive systems. 
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1. Introduction 

Transforming a real- world system into a deterministic model (e.g., a set of 
dynamic equations) often leads to errors whose effect may greatly influence 
the process and can not be neglected. 

If the nature of errors is random or probabilistic, then it is possible to 
adopt a stochastic framework; however if the underlying structure is not 
probabilistic, for instance due to subjective modeling choices, a different 
formalism is required [8]. 

Fuzzy theory, first introduced in 1965 by L.A. Zadeh [10] . appears the 
most natural choice to model complex systems affected by vagueness. 

In the literature many approaches have been introduced: in [T2| H3] fuzzy 
numbers are used to model uncertain quantities (e.g., different beliefs or 
opinions [20]), considering also the fuzzy extension of traditional arithmetic 
operations; a different approach is to synthesize controllers based on fuzzy 
rules [131 IT5l [16l [T7] in order to control complex real- world systems encoding 
the experience of human operators (e.g., a parking system for a wheeled 
mobile robot [18] or the navigation of a mobile robot [19]); a different, more 
formal approach is to consider a dynamic fuzzy equation or system, both in 
the continuous [2D El [221 E3J [21] and discrete-time [SJ fashion. 
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The role of Fuzzy theory appears of particular practical utility when hu- 
mans are directed involved; for instance in [TTj . systems with crisp dynamic 
and fuzzy state have been studied, and the framework has been applied to 
critical infrastructure related problems, where the only information available 
is that provided by operators and stakeholders and is therefore linguistic and 
vague. 

In this paper the stability of discrete-time linear and stationary fuzzy 
systems where both the coefficients and the state are described by means 
of fuzzy variables is studied, allowing to represent uncertainty to a greater 
extent with respect to the approach in [TT] . 

In the continuous time case, the approach of representing the fuzzy system 
as a family of Differential Inclusions [H El IE] has been successfully adopted in 
[HE1IZ]; in this paper we will follow a similar path, introducing the framework 
of Fuzzy Difference Inclusions (FDI). 

The paper is organized as follows: after some preliminary definitions, 
Section 2 introduces the Discrete-time Fuzzy Systems and Section 3 specifies 
them as a family FDIs, while Section 4 address the stability of fuzzy systems 
described as a family of FDIs; Section 5 further specifies the results for 
linear systems, while Section 6 characterizes the solution set of positive linear 
systems; finally some conclusive remarks are collected in Section 7. 

1.1. Preliminaries 

In the following, vectors and vectorial functions will be represented by 
boldface letters while scalars and functions with scalar codomain will be 
represented by plain letters. Moreover, to avoid confusions, x will denote a 
vector with fuzzy entries, while crisp (i.e., non-fuzzy) vectors will be denoted 
by z. 

Let I p denote the p x p identity matrix and let c p be a vector of p com- 
ponents, each equal to c. 

Let 1R, N be the set of reals and integers, respectively and IR + ,N + be the 
set of nonnegative real and integer numbers, respectively. Let be the 
space of nonempty compact convex subsets of R . Let be the open unit 
ball in R N . 

Given a space X and a particular distance d defined over such a space, 
the Hausdorff separation and the Hausdorff metric for two sets A,5cX are 
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given by: 



Pd X (A, B) = sup{d(a, b) : a G A,b G B} 
Pd>x (A, B) = max{p* dtlL (A, B),p* d<x (B, A)} 



(1) 
(2) 



In this work it will be considered the following distance in R 



JV 

d RN (zi,z 2 ) = } j d M .(z 1 i,z 2 i) 



(3) 



i=i 



where d^(-, •) is the euclidean distance in R and z ly z 2 G M. N . 

If a matrix M is non-negative (non-positive), i.e., it has only non-negative 
(non-positive) entries, write M > (M < 0). If B — A is a nonnegative 
matrix, write B > A. 

2. Fuzzy Discrete-time Systems and Inclusions 

The key idea of fuzzy theory is to extend traditional set theory by allowing 
an element to partially belong to a set. Therefore, a fuzzy subset of R is 
defined by a membership function x : R — > [0, 1] which assigns to each point 
p G R a grade of membership in the fuzzy set [10J. Notice that, when dealing 
with fuzzy sets, it is very common to denote with the same symbol the fuzzy 
sets itself and the membership function that defines the set, because the 
membership function univocally determines the set and vice-versa. 

Let the p-membership x(p) be defined as the grade of membership of p G M. 
in the set x. In the following, the value assumed by a given fuzzy variable at 
time instant k will be denoted by x(k); in this case the p-membership of the 
set x(k) is denoted by x(p, k). 

For each a G (0,1], the a-level set [x] a of a fuzzy set is the subset of 
points p G K with membership grade x(p) > a. The support [x}° of a fuzzy 
set is defined as the closure of the union of all its a-level sets (see Figure [T|. 

The use of a-levels, as shown in Figure [TJ allows to treat fuzzy sets as a 
set of nested real intervals; in the following the evolution of a fuzzy dynamic 
system will be evaluated level-wise by considering for each a, the evolution 
of a system whose state is described by an interval (i.e., the corresponding 
a-level). In [8] E is introduced as the space of all fuzzy subsets x of R such 
that: 

1. x maps R onto [0, 1]; 
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Figure 1: Different a-levels of a triangular-shaped fuzzy membership function. 



2. [x]° is a bounded subset of R; 

3. [x] a is a compact subset of R for all a G (0, 1]; 

4. £ is /uzzy convex, that is: a;(0p + (1 — </>)g) > min[x(p),x(q)] for all 
p, q e R 

These condition are equivalent to require that the membership functions 
are compact, that their shape is composed of a monotone nondecreasing and 
a monotone non increasing part (e.g., a triangle), and the a- levels are nested, 
i.e., those with greater values of a are contained into those with smaller a 
(all the intervals are contained in the support). With this formulation the 
fuzzy sets in E are often called Fuzzy Numbers (FN) [T2l [TB"] . 

A triangular fuzzy number (TFN) x G E, in particular, is described by an 
ordered triple {xi, x c , x r } G R 3 with xi < x c < x r and such that [x]° = [xi, x r ) 
and [x] 1 = {x c }, while in general the a-level set is given, for any a G [0, 1] 
by: 

[x] a = [x c - (1 - a)(x c - xi), x c +(l- a)(x r - x c )] (4) 

An illustrative example of triangular fuzzy numbers is reported in Figure 
[2j the figure shows two triangular fuzzy numbers with different levels of 
ambiguity, representing for example the codification of the statement "an 
high wall" . Specifically, the blue internal triangle represents the fuzzy number 
corresponding to the measure of 4 meters with a given degree of ambiguity 
(i.e., the width of the base of the triangle); the red external triangle still 
represents the same fuzzy number 4, but with a greater ambiguity. 
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The triangular representation is not the sole available alternative; as de- 
picted in Figure [3] many other shapes are possible, and more complex is the 
shape, more descriptive is the resulting fuzzy number (i.e., the vagueness is 
better characterized). For instance the existence of a plateau for a given in- 
terval represents complete indeterminacy for that interval, or an asymmetry 
with respect to the peak may represents different beliefs in the entity of the 
values that are smaller or bigger than the "nominal" value (e.g., the best and 
worst cases). However, in many applications, the Triangular Fuzzy Numbers 
are the most used, because they can be described by the triple of the abscis- 
sae of their vertices ({2,4,6} and {3,4,5} in the case of the Triangular FN 
in Figure BJ); moreover, as shown in Figure [TJ a crisp value is obtained for 
a — 1, hence the level with the strongest belief collapses into a point. 

More complex shapes, as shown in Figure |3j allow to better characterize 
the uncertainty: the leftmost Fuzzy Number represents a situation where un- 
certainty rapidly decreases while approaching the peak value; the rightmost 
Fuzzy Number, due to its trapezoidal shape, models the case where a single 
value with maximum belief can not be found. Notice further that, as stated 
before, the shape of a FN needs not to be symmetric (see the central TFN 
in Figure [3]), thus allowing to represents different beliefs on the left and right 
spread of uncertainty with respect to the value associated with the maximum 
belief. 

1 - 
0.9 - 
0.8 - 
0.7 - 
0.6 - 
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- 
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Figure 2: Two triangular fuzzy numbers representing the same value "about 4" , although 
with different uncertainty. 
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Figure 3: Examples of fuzzy numbers in the space E. 



In order to measure the "distance" between two fuzzy sets, the space E 
is equipped with the metric [8]: 

dE{xi,x 2 ) = sup{|a:i(» - x 2 (p)\ : peR} (5) 

which measures the largest difference in the membership grades of two fuzzy 
sets x\ and x 2 ; clearly d^(xi,X2) G [0,1]- Figures [1] and [5] provide examples 
of computation of such a distance; it is worth to notice that, unless the 
maximum a-level (i.e., a = 1) of one of the two FNs coincides with an interval 
where at least a point in the other FN is non-zero, then ^(^1,^2) = 1 (see 
Figure [4]); when such a condition is verified, the metric assumes a value that 
is inversely proportional to the degree of overlapping of the two FNs (see 
Figure |5]); finally, notice that the distance becomes zero if and only if the 
shape of the two FNs coincide. 

The following equation correlates the distance with the a-sets [S] , consid- 
ering the euclidean distance d in R: 

d E (x 1 ,x 2 ) = sup{p d)R ([a;i] a , [x 2 ] a )} x h x 2 G E (6) 

In the following the concept of level-wise convergence of a sequence of 
fuzzy numbers will be defined, that is, for a fixed a G [0, 1], the convergence 
of the a- levels of the FNs in the sequence. 

Let {x n } be a sequence on E, then {x n } converges level- wise to x G E if, 
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for all a G (0, 1]: 

pd R ,B,([ x n] a , M a ) -> as n oo (7) 

Define 

$ = {iGE: + (1 - 4>)q) > 4>x(p) + (1 - (8) 

for each p, q G G [0, 1]. 

In order to consider vectors of N components, each being a FN, the space 
K N has to be characterized; to this end K N is equipped with the following 
metric: 

N 

d E iv(x,y) = ^2d E (xi,yi) (9) 

i=l 

where x = . . . , x n ] T and y = . . . , y n ] T , x, y G E n , while Wi = 1,...,N 
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Figure 5: Example of distance calculation for two TFNs x\ — {2, 3, 4} and X2 — {0, 3, 8}. 
The distance d%(xi, x^) is equal to 0.8. 



x h yi E E. 

Define the a-level of a vector of FNs x e E^ as the set of vectors z e M. N 
such that, Vi = 1, . . . , N z% belongs to the a-level of i-th component Xj. 

It is worth to notice that many applications based on fuzzy theory often 
handle uncertain models implicitly, by considering approximations of the 
algebraic operations extended to fuzzy numbers; however these approaches 
are often limited to triangular shapes, due to computational constraints and 
due to the complexity of the approximations if the shape is not triangular. 
Other approaches are typically aimed to synthesize controllers based on fuzzy 
rules p3l IT5l [16l [T7] in order to control complex real- world systems encoding 
the experience of human operators (e.g., a parking system for a wheeled 
mobile robot [18] or the navigation of a mobile robot [E]). Within all these 
approaches, however, it is not possible (or very hard) to obtain a closed form, 
neither in the simplest cases, and it is therefore very difficult to study the 
stability of the fuzzy systems. 

In order to clarify the limits of simulation-based techniques, consider for 
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instance the scalar system 



x(k + 1) = mx(k) + q (10) 

where x(0) = {x /, x 0c , x 0r }, c = {q,c c , c r } and q = {qi,q c ,Q r } are triangular 
fuzzy numbers. In the simulative perspective, the extension of the sum to 
TFNs is immediate; indeed, such an operation is linear and the resulting 
fuzzy number is also triangular [12]. The sum of two TFNs = [a h a c , a r ] and 
b = [bi,b c , b r ) is given by 

a + b = [ai + b[, a c + b c , a r + b r ] (11) 

Analogously, the multiplication by a real scalar /3 is given by 

/3a — [Pax, /3a 2 , /3a 3 ] (12) 

Triangular fuzzy numbers, therefore, are closed with respect to sum and 
product with a real scalar [12J. 

When addressing the product of two TFNs, however, there is the need to 
introduce approximations, since the product of two TFNs is in general not 
triangular; in the literature many approximations of the product have been 
given, and the most simple, in the case where a, b > (i.e., a/, b\ > 0), is 
given by: 

a* b = [ai * bi, a c * b c , a r * b r ] (13) 

However, in this way, besides loosing expressivity (because only triangular 
numbers can be typically handled by simulators), it is not possible to analyze 
and formally characterize the stability of the system. 

Moreover the approximated methods fail to address more complex situa- 
tions, where, besides the products, a general nonlinear fuzzy-valued function 
with fuzzy argument has to be computed. 

In the following we will consider Complete Discrete-Time Fuzzy System 
(C-DFS), defined as follows: 

x(A; + l) = H(x(A;),A;), x(0) = x (14) 



where H : N+ x E N -»■ E N and x,x G E^. Notice that for System Q, 
even if the initial conditions were crisp, due to the fuzziness of the dynamic 
the state of the system is composed of fuzzy variables. 
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However, as noted before, the dynamic of system (14) generates evolutions 
that are in general complex and not triangular, and the level-wise evaluation 
of such a system is particularly hard, as will be explained in the following; 



this makes unfeasible any direct calculation of (14). 

An elegant solution can be achieved in the framework of Difference In- 
clusions E]. 

3. Fuzzy Difference Inclusions 

In the continuous time case, the approach of representing a fuzzy system 
as a Fuzzy Differential Inclusion (FDI) that is, a family of difference inclu- 
sions [TJ [3] (one for each a-level) , has been successfully adopted in [U [21 E] • 



Defining [H(x(fc), k)} a = Q([x(k)] a , k, a), the C-DFS (Q can be rewrit- 
ten as the FDI: 

x a (fc + 1) £ Q([x(/c)] a , k, a); x Q (0) £ [x ] a ; a£j (15) 

where J = [0,1], Q is an open subset of M N x N + containing [xo] a x {0} and 
Q : Q x J — y is a continuous and set valued map. This latter condition is 
equivalent to require that the set-valued function Q([x(/c)] a , k, a) is defined 
for a neighborhood of the initial conditions, and its image is compact and 



convex. Notice that in x a (k + l),x Q (0) £ R N and hence [x(jfe + l)] a is 
the set obtained considering all the attainable x. a (k + 1), starting from all 
the x. a (k) £ [x(A;)] a . 

From now on, it is assumed that all maps are proper, that is have nonempty 
images of points in their domain. 

The idea of considering difference inclusions becomes clearer while con- 
sidering a simple example; consider the system 

x(k + 1) = mx{k) x(0) = x (16) 

where m = {mi, m c , m r } and xq = {xoi, x$ c , xq t } are fuzzy. As stated before, 
the product does not yield to a triangular fuzzy number and is very hard 
to evaluate; this is even more true in a level-wise representation. Indeed, 
for each a-level we have to take into account the product of two intervals, 
[rn a ,m a ] and [x a (k),x a (k)} and the approach used for simple DFS system 
cannot be adopted. A solution is to consider, for each a-level and for each 
point belonging to the interval [m a ,m a ], a level-wise representation of a 
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standard DFS, thus obtaining the difference inclusion 



x a (k + 1) G [m] a [x(k)} a (17) 

In order to proceed there is the need to recall some basic notions on 
difference inclusions and their stability. 

3.1. Difference Inclusions 

A Difference Inclusion [5], [6] is given by 

z(k + 1) G P(z(&), k), z(0)Gz (18) 

where z C R , f2 is an open subset of R x N + containing z(0) x {0} and 
P : O — > is a continuous set valued map. Notice that the dynamic P is 
set valued and also the initial condition is a set. 

Notice further that, in a very general perspective, difference and differen- 
tial inclusions are often required to be upper or lower semicontinuous, since 
it is possible to study the behavior of inclusions with discontinuous right-end 
side; however such an extension is out of the scope of this work, and the 
reader may refer to [HI El HI E] for a more general approach. 

Let Z k (R N ) be the set of continuous functions / : R N x {0, . . . , k} -> R N . 
We need to define the following sets: 



Set of Solutions: the set of solutions S(z ,A;) of the inclusion (18) 



on {0, . . . , k}, is the set of all solutions of the inclusion (18) from time 
step until k; clearly S(z , k) E Z k (R N ). 

Attainable Set: the attainable set, for a given step k, is the set of 



values that the solutions of the inclusion ( 18 ) assume in k and is defined 

as 

A(z 0l fc) = {z(fc):z()GS(z ,fc)} (19) 



A set M is stable for the Inclusion (18) if for all e > there exists 5 = 
5(e) > such that z G M + 5M N implies z(Jfe) G M + eB^ for all k G N + and 
for every solution in §(z , k). 

The above definition may be rephrased as zq G M + 518^ implies that 
Pd R2V ^(A(z ,A;),M) < e for all k G N+, where p dRjV)R iv(-, •) is the Hausdorff 
distance in R-^. 

If p disJV jR jv(A(z , k), M) — > as k — > oo and M is stable then it is asymp- 
totically stable. 
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4. Stability of Fuzzy Difference Inclusions 



Let us now discuss the characteristics of a C-DFS ( 14 ) expressed as a FDI 



(15). 



Denote the set of solutions of an inclusion in the family ([15j) (i.e., for 
a given a G [0,1]) by S a (xo,^) and the attainable set by A a (xo,/c), while 
§^(x , k) and A-^(x , k) are the set of solutions and the attainable set for the 



whole family (15). 



Note that, during the evolution of a FDI, it is not in general verified that 
the convexity of the a-levels is preserved; therefore the stability of FDIs is 
often studied in the space of fuzzy sets with not necessarily convex level 
sets [U I3T] ;however in this work we relax the semicontinuity assumption, 
supposing that the dynamic is continuous, hence avoiding the issue of non- 
convex level sets [3T] • 

The following Stacking Theorem [8], whose statement is reported below, 
characterizes the elements of D, but can be easily generalized to the vectorial 
case or to any Banach space. 

Theorem 4.1. Let {Yp C M> : < (3 < 1} be a family of compact subsets 
satisfying 

1- Y p C Y a for < a < (3 < 1 

2. Yg = H^Yp. for any nondecreasing sequence fa —> (3 in [0, 1] 

Then there is a fuzzy set «6D such that [u] a = Y a . If the Yp are also convex, 
then u G E. Conversely the level set of any mgE are compact and convex 
and satisfy these conditions. 

The above theorem states that a family of compact and closed nested 
intervals coincides with a fuzzy set in D, and the property can be trivially 
extended to the vectorial case (if the property is true for each state variable); 
moreover a fuzzy set in E^ is always composed of level sets that are non 
increasing in a. 

There is the need to introduce the following stability definitions, that 
extend to the fuzzy domain the definitions provided for difference inclusions. 
Notice that, although xo is a point in E^, the stability definitions involve 
the fuzzy attainable set, therefore the stability is defined for a set U C E^. 



A set U C E, N is stable for System (15) if for all e > there exists a 5(e) 
such that 

x G U + -» A / (x ,fc) G U + eB^ (20) 
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for all k G N + , where is the open unit ball in E . 
In other terms 

P* eJViE ^(x ,U) < 5 p dEN<EN (Af(x ,k),l]) < e (21) 

for all k G N + , where p* d nE n(-,-) and Pd EN ,E N (^ ') are the Hausdorff distance 
and the Hausdorff separation in E^, respectively [8], based on a distance 
d^N, defined in 

If p^ N ^N (AjJxq, k) , U) — > as k — > oo and U is stable, then it is asymp- 
totically stable. 

The following theorem characterizes the structure of the set of solutions 



and the attainable set of a FDI (15). 



Theorem 4.2. Let UcD" and suppose that [U]° C R N is bounded. The 



attainable sets A Q (xo,/c) of the family of inclusions (15) are the level sets of 
a fuzzy set A^(x , k) and the sets of solutions S a (x , k) are the level sets of 
a fuzzy set §^(x , k). 

Proof 4.3. The set§ Q (xo, k) is composed of solutions x a (•) such thatx. a (k') G 
A a (xo, k') for each k' G {0, . . . , k}. 

To show that the sets §«(xo, k) for each a G [0, 1] are compact we need to 
show that A Q (x , k') is compact for each k! G {0, . . . , k}. It is a well known 
result that compactness is invariant to continuous transformations (see ISO] 
for instance); since [x(0)] Q is compact, it follows that A Q (x ,l) is compact. 
Iterating it is verified that all the A a (xo, k') are compact and hence S a (xo, k) 
is compact. 

Notice that Q([x(A;)] a , k, a) C Q([x(A:)] , k, 0), [x] Q C [x]° for all a G 
[0, 1], therefore the sets § a (x , k) are nonincreasing in a. 

Let {Pi} be a nondecreasing sequence in [0,1] converging to (5 G [0,1] 
and consider the sequence {S^xo, k)}; clearly the intersection of the first j 
elements in the sequence coincides with S y g i (xo, k). Therefore the intersection 
of all the elements in the sequence gives 

U§ ft (x o ,A;)=Mxo,£0 (22) 

i 

where, for infinite sequences, the equivalence holds in the limit of the infinite 
intersection. A similar argument proves the result for the attainable sets 
A a (x , k). 
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It is now possible to apply the Stacking Theorem (4-1) on the space Zfc(M 7V ) 
for S J g(xo, k) and on the space V> N for A a (xo, k), proving the statement. 



The following theorem characterizes the stability of a FDI (15). 

Theorem 4.4. Let U C Ib N and suppose that [U]° C M. N is bounded and the 
inclusion 

[x]°(* + 1) G Q([x]°(fc), k, 0); [x]°(0) = [x ]° (23) 
is (asymptotically stable); then the FDI |73ty (asymptotically) stable. 



Proof 4.5. From Theorem 4-2, the FDI has a fuzzy attainable set A^(x , k) 
which coincides with [A^(x ,k)]° , because the other level sets are nested in 
the support. 

Since the difference inclusion (23) is stable (there exists a stable set M), 
then A^(xo, k) is bounded and the family (15) is stable. 
If the inclusion (23) is asymptotically stable, then 

p dRjVjM iv(Ao(x , k), M) ->■ as k -> oo (24) 

Therefore, since [A* (xo, k)] a C A^(x , k) for < a < 1 it follows that 

p d ^N([A f (x , k)] a } M) -> as jfe ->■ oo (25) 



and tae /uzz?/ system (15) is asymptotically stable. 

The above theorem proves that the stability of the level-wise representa- 
tion of the support of the fuzzy system (i.e., a = 0) is a sufficient condition 
for the stability of the whole fuzzy system. 

We can now provide some further characterization in the case of linear 
and stationary FDIs. 



5. Linear and Stationary Fuzzy Difference Inclusions 

In this section we will consider Linear and Stationary C-DFS, defined as 
follows: 

x(Jfe + 1) = Hx(Jfe), x(0) = x (26) 

where H is a n x n fuzzy valued matrix, i.e., hij G E is a fuzzy number and 
x, x G E^. Let [H] a = [H tt , H„] for any a G [0, 1] be an interval matrix, 
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where H a , H Q denote the matrices whose elements are, respectively, the lower 
and the upper end points of the interval. 



The linear and stationary C-DFS (26) can be rewritten as the FDI: 
x a (k + 1) G [H] a [x(k)] a ; x Q (0) G [x ] a ; < a < 1 



(27) 



It is possible to further specify the results of Theorem 4.4 



Corollary 5.1. Let a Linear and Stationary FDI (27); then: 



1. ifVU G [H]° is stable then the fuzzy system (21) is stable. 

2. if WU G [H]° is asymptotically stable then the fuzzy system (27) is 
stable. 

Proof 5.2. If the first condition is verified, the attainable set [A' (xo, k)]° is 
bounded and the inclusion (23) is stable. From Theorem 4-4- we have that 
also the linear and stationary FDI (21) is stable. For asymptotical stability 
note that if the second condition is verified, then the fuzzy System (21) is 



asymptotically stable, again by Theorem 4-4 



The above corollary, therefore, requires the inspection of all the dynamic 
matrices in the interval [H]° to prove the stability; under rather general hy- 
potheses, the following result provides more operative necessary conditions. 



Corollary 5.3. Let a Linear and Stationary FDI (21) and suppose that the 
crisp matrix H > and that, for all i = 1, . . . , N 



N 

E ( H °) 



< 1- (H 



0)ii 



(28) 



then the FDI (21) is asymptotically stable. Analogously the result holds if 
H < and 

N 

j=l;j¥=i 



> 



:h ), 



(29) 



Proof 5.4. Since H > it follows that for each matrix U G [H]q, U > 0; 
moreover H > is asymptotically stable, since by Gershgorin circle Theorem 
19$, its eigenvalues lie, in the complex plane, in the union of circles centered 
in G [0, 1] with radius equal to r 1 ^ = J2f=i-j^i 



H 



0)ij 



< 1 



H 



0)ii- 
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Since H > U, for all U G [H]°, we have that 

max(uii + r\j) < max((H )« + rL) < 1 (30) 

i i u 

and each matrix U is asymptotically stable in the discrete-time sense. Ac- 



cording to Corollary 5.1 the fuzzy System (21) is stable. The proof for the 



negative case is analogous. 

Notice that the condition H > includes several systems of practical interest 
as, for example, the class of positive system whose Gershgorin circles are 
contained inside the unit ball. These systems find large applications in vague 
data applications (see, for instance, [2T| 128])- 

The following corollary applies to a generic linear and stationary FDI. 



Corollary 5.5. Let a Linear and Stationary FDI (27), and consider the 
interval matrix [H]° = [H ,Ho]. 

Define H 0c = |(H + H ) and AH = ^(H - H ). Let [r)° = [r ,f ] and 

= [i , io] be defined as follows: 



r = min {x T H 0c x — |x| T AH|x|} (31) 

IN|2 = 1 

f = max {x T H 0c x — |x| T AH|x|} (32) 

IM|2 = 1 

i = min {xj (H 0c — Hq c )x 2 — AH o \x1X2 — x 2 x\ |} (33) 

||tasi,a:a]||2=l 

i — min {x^(H 0c + Hq c )x 2 - AH o \xix% - x 2 x\ |} (34) 

||[xi,a!a]||2=l 

where AoB = a>ijbij is the scalar product of the matrices A and B and 

[xx, X2\ is the stack vector composed of x\ and x 2 (where the transposition of 
vectors has been omitted for brevity). 
If the following condition holds true 

max{||r + ji \\, ||r + ji \\, \\f Q + ji IUI f o + jio\\} < 1 (35) 



where j is the imaginary unit, then the FDI (21) is asymptotically stable 



Proof 5.6. In [25] it is proven that, given a square interval matrix [H , H ] 
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for each matrix H Qi G [H_ ,H ] the eigenvalues Xq^ ofH 0i are such that: 



r < Re[X 0ij } < r 
io < Im[X 0ij ] < i 



(36) 
(37) 



Therefore, if condition (35) holds true each matrix H 0i is asymptotically 



stable in the discrete-time sense. Hence, according to Corollary 5.1 the fuzzy 



System (21) is stable. 



Notice that the above operative conditions are just sufficient conditions, 
and that such conditions do not address the case where the system is stable 
but not asymptotically stable, because these methods can not determine the 
multiplicity of the eigenvalues on the unit circle of the matrices in the interval 
[Hf. 

The following corollary provides a simple sufficient condition to determine 
the simple stability of the fuzzy system. 



Corollary 5.7. Let a Linear and Stationary FDI (21) and suppose that there 
exists a transform matrix T such that 



T _1 H n T 



or 



T^HqT 



H? 



(38) 



If the conditions of Corollary \5.3 hold, but inequality (28) is not strictly 
satisfied (i.e., ^. =1 -^(Hq)^ < 1 — (H-o)u) and inequality (28) holds 
strictly for the N — 1 x iV — 1 matrix H , then the FDI (21) is stable. 



If the conditions of Corollary 5.3 hold, but inequality (29) is not strictly 



satisfied and condition (38) holds for H Q (inequality (29) holds strictly 



for matrix Hq^, then the FDI (21) is stable. 



If the conditions of Corollary 5.5 hold for the interval matrix [Hq, Hq] 



and condition (38) holds for both H and H 0; then the FDI (21) is 
stable. 

Proof 5.8. From Gershgorin circle Theorem fS/ the eigenvalues of each ma- 
trix U G [H]° are such that ||Aj|| < 1 and Re[Xi] > 0. 

In order to prove that the system is stable (but not asymptotically stable), 
we need to prove that the multiplicity of eigenvalues with ||Aj|| = 1 is exactly 
1. 
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Clearly, under the first set of hypotheses matrix Corollary 5.3 holds strictly 
for the reduced matrix Hq, and therefore it has eigenvalues within the open 
unitary circle in the complex plane; due to the triangularity of the trans- 
formed matrix, the statement is verified. An analogous argument proves the 
statement for the other cases. 

The following section describes the practical representation of a linear and 
stationary Discrete-time Fuzzy System, and is the discrete-time analogous of 
the results in [8j [2] . 

6. Representation of Linear Fuzzy Difference Inclusions 

In order to study the practical representation of a linear and stationary 
Discrete-time Fuzzy System there is the need to discuss the solution set of 



System (27). Consider a linear crisp system z(k + 1) = G(/c)z(fc),z(0) = z . 



The solution of such system is given by 

z(k) = ®(k)z (39) 
where $(/c) satisfies the matrix difference equation 

$(k + 1) = G(fc)$(fc); $(0) = Jtv (40) 
Clearly, if the system is also stationary, $(&) = G k . 



Analogously to the crisp case, consider the fuzzy version of System (39); 



then, equation (|40j) can be interpreted as the family of difference inclusions 
$ Q (fc + l) G [H(k)] a <5> a (k); $ a (0)=I N ; 0<a<l (41) 

Denote 

$ a (h) = {Y(k) : Y{k + 1) = U{k)Y{k), Y(0) = I, [/(-) G [#]«(•)} (42) 



From the Stacking Theorem 4.1, the are the level sets of the D iVxAr 

valued fuzzy function $(■). If the matrix H is stationary and [H] a is an 
interval matrix for each a G [0, 1], then $(•) is E iVx7V valued, since an interval 
matrix is a convex set in M^^. In this case U(k) = U and Y(k) = U k , hence 

$ Q (Jfe) = {U k :U E [H} a } = [H k } a (43) 
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In this case denote the fuzzy set as H k . 

In the following theorem it is proved that, if is nonnegative, then, for 
each level set, the evaluation of bounds of the set of solutions is simplified. 

Notice that, within the proposed approach, the solution is not a single 
trajectory, but it is defined as a set; there is the need to consider, therefore, 
the attainable set and the set of solutions. 

In order to derive a framework with a real applicability, the following 
theorem proves that, limiting the scope to the class of linear and stationary 
FDI with non-negative entries, the evaluation of the upper and lower bounds 
of the set of the solutions is extremely simplified. 



Theorem 6.1. Let a linear and stationary FDI (21) and suppose that the 



crisp matrix H° > 0; then the following holds for each a E [0, 1] 

§ a (x , k) = (H a ) k x a0 ; § a (x , k) = (H a ) k x a0 (44) 

where x a0 and x a0 are the left and right bounds of [xq]", respectively and 
§ Q (x ,A;), §> a (x ,k) are the left and right bounds of the set of solutions 



§q,(x , k) of the single Inclusion in the family (15). 



Proof 6.2. We have to show that, given two matrices R and G, with G > 
R > such that 

X(k + 1) = RX(k), X(0) = I 
Y(k + 1) = GY(k), Y{0) = I 

it follows that Y(k) > X(k),k > 0. Clearly, X(k) = R k and Y(k) = G k . 
Since G > R > the statement is proved. Note that if H > then the level 
set matrices and H k are positive matrices for each a G [0, 1], and when 
multiplying interval vectors the order of endpoints is preserved. We have 
therefore that Sq,(x , k) = [§ Q (x , k), S a (x , k)] and the theorem is proved. 

The theorem proves that, under the hypotheses, for each a-level the upper 
and lower bound of the set of solutions can be calculated independently as if 
they were two crisp systems with dynamic matrices if a and H a , respectively; 
therefore, for practical applications, it is sufficient to consider the bounds to 
adequately represent the solution set. 
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7. Conclusions 



In this paper the stability and representation of linear and stationary 
fuzzy systems where both the dynamic coefficients and the state variables are 
described by means of fuzzy numbers is addressed, providing some stability 
conditions. 

Further works will be devoted to find stability conditions for time varying 
and nonlinear fuzzy systems. 
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